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Lau, Berciu and Sawatzky reply: In the preceding 
Comment , Lee and Lee bring to attention their inter- 
esting variational calculation for the single band Hubbard 
model further reduced to a t-t'-t"-J model 0. Its main 
result was to reveal new low-energy one-hole states called 
spin-bags (SB), possibly forming a continuum. SBs con- 
sist of a quasiparticle (QP) plus a spin- wave excited in the 
AFM background, and are found to cross below the QP 
band in some regions of the Brillouin zone (BZ). Based 
on this, Lee and Lee claim that the SBs explain, within 
a one-band model, the spin- |_ polaron that we recently 
found in a three-band model [3[ . They conclude that our 
claim that this model reveals physics that cannot be de- 
scribed within one-band models, i.e. in the framework of 
Zhang-Rice singlets (ZRS), is not justified. 

While superficial similarities exist between the SB and 
the spin- 1 polaron, we disagree that they describe the 
same physics, on several grounds: 

(i) In the Supplemental Material of Ref. , we ruled 
out the possibility that the spin-| polaron is a spin-i 
polaron plus a free magnon, because its band lies below 
the continuum describing such states. It can be roughly 
thought of as a bound-state of a spin-i polaron and a 
magnon, with a very distinct local spin structure around 
the charge. The existence of such bound states, which 
might be a better analog of our spin-| polaron, is not 
analyzed for the one-band model, in Refs. 0,0; 

(ii) As shown in our Fig. 2, the spin— | polaron's 
band has significant dispersion, comparable to that of 
the spin— i polaron 0. In contrast, the low-energy edge 
of the SB continuum is rather flat throughout the BZ, 
see Fig. 1(b) of Ref. 0. This striking difference in their 
spectra is likely an indication of a very different nature of 
the two types of low-energy states. There is currently no 
evidence that the two models have comparable dispersion 
for spin excitations, regardless of their nature. 


A second point raised in the Comment is that if a 
ZR-like state is built from a superposition of configu- 
rations like that of Fig. 3a, AFM correlations on the e 
and d bonds are similar to those calculated in Ref. 0. 
This is taken as proof that the spin-i polaron is simi- 
lar to the ZR-based QP state, as well. First, Fig. 3a 
is for a state of momentum ( § , § ) , so naive ^ rotations 
lead to a state with an ill defined momentum. In fact, 
even though these bonds are related by the exact P x ^ y 
symmetry of our Hamiltonian, the quoted values are not 
equal; this is wrong. In any case, the fact that bonds 
rather far from the hole show robust AFM correlations 
is hardly surprising. The key observation in our model is 
the strong FM correlation between the spins neighboring 
the hole, which points to the three-spin polaron (3SP) as 
the proper framework to understand the spin-i polaron 
and the inner core of the spin-| polaron. Since the 3SP 
can be written as the sum of singlets between the hole 
and each of its neighboring spins [3j , it does have a finite 
overlap with a ZR state Q . Its additional degrees of free- 
dom, however, allow it to describe correlations beyond 
those possible in a ZR-based model. This invalidates the 
Comment's claim that a low-energy non-bonding state is 
the only signature of breakdown for a one-band model. It 
is very important, in this context, to also point out that 
the model used in Ref. is a further simplification of 
the ZR scenario - the O states are no longer present and 
a discussion of the spin correlations around an O hole 
becomes meaningless. 

For all these reason we remain convinced that both 
the spin-i and spin-| polarons in our 3 band model are 
quite different objects from the QP and SB obtained from 
a single band description However, caution is necessary 
since comparisons between a variational solution based 
on mean-field and our exact diagonalization (ED) for a 
finite cluster may be misleading. If ED results for a one- 


(iii) While the spin-| polaron band crosses below the band model revealed similar low-energy spin-| states, 


spin-i polaron band in certain regions, just like the SB 
continuum is below the QP band in certain regions of 
the BZ, a careful comparison shows yet more differences. 
In our model, this happens in two separate regions, cen- 
tered at (0,0) and (7r,7r). In the variational solution for 
the one-band model, this happens in one larger region 
centered at k = (ir, it) which, coincidcntally, is the AFM 
order vector. The difference is most clearly visible along 
the (ir, 0) — (0, it) cut, where we find no crossing whereas 
the variational calculation predicts the QP as the low- 
energy state only near (f , f )• If bound-states were found 
in the one-band model, the comparison would be worse 
since this would further increase the crossing region. 

Such differences result in very different physics, eg. at 
the nodal point. While the vanishing quasiparticle weight 
at (0, 7r) is explained as being due to the SB state in 
the one-band model, we find Z = here because of the 
orthogonal reflection parity between the lowest electron- 
removal state and the lowest spin-i eigenstate 0. 


and FM correlations between the spins sandwiching the 
hole, our position would have to be reconsidered. To our 
knowledge, the former is not the case and the latter is 
not possible for a one-band model. 
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